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ABSTRACT 

^-5 ; 

We investigate the physical conditions of the interstellar medium and stellar 

components in the regions of the southern Galactic star forming complexes as- 

^; sociated with IRAS 10049-5657 and IRAS 10031-5632. These regions have been 

I ! mapped simultaneously in two far infrared bands (Ag// ~ 150 & 210 /im), with 

^ ! ~ 1' angular resolution using the TIFR 1-m balloon borne telescope. Spatial dis- 

c/3 , tribution of the temperature of cool dust and optical depth at 200 /im have been 

obtained taking advantage of the similar beams in the two bands. The HIRES 

processed IRAS maps at 12, 25, 60 and 100 /im have been used for comparison. 

^ ! Using the 2MASS near infrared sources, we find the stellar populations of the em- 

■rj- ! bedded young clusters. A rich cluster of OB stars is seen in the IRAS 10049-5657 






region. The fits to the stellar density radial profile of the cluster associated with 
IRAS 10049-5657 has been explored with inverse radius profile as well as the 
O I King's profile; the cluster radius is obtained to be ~ 2 pc. The source in the 

Q i cluster closest to the IRAS peak is IRA-7 which lies above the zero age main 

sequence curve of spectral type 05 in the colour-magnitude diagram. Unlike 
k> [ IRAS 10049-5657, a small cluster comprising of a few deeply embedded sources 

H ' is seen at the location of IRAS 10031-5632. Self consistent radiative transfer 

modelling aimed at extracting important physical and geometrical details of the 
two IRAS sources show that the best fit models are in good agreement with the 
observed spectral energy distributions. The geometric details of the associated 
cloud and optical depths (rioo) have been estimated. A uniform density distri- 
bution of dust and gas is implied for both the sources. In addition, the infrared 
ionic fine-structure line emission from gas has been modelled for both these re- 
gions and compared with data from IRAS-Low Resolution Spectrometer. For 
IRAS 10049-5657, the observed and modelled luminosities for most lines agree 
to within a factor of four while for IRAS 10031-5632, we find a discrepancy of a 
factor of 100 and it is likely that some basic assumptions of the model are not 
valid in this case. 
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Subject headings: stars: formation - infrared: ISM - ISM: H II regions - ISM: 
individual (IRAS 10049-5657, IRAS 10031-5632) 



1. Introduction 

HII regions are astronomical sources that represent early stages of deeply embedded 
high mass (O or early B) stars. Their study can provide vital information about high 
mass star formation as well as their interaction with the parent molecular cloud. Being 
deeply embedded in interstellar cloud including the dust component, almost all of their 
energy is absorbed and re-emitted in the infrared wavebands. IRAS 10049-5657 (G282.0- 
1.2) and IRAS 10031-5632 (G281. 6-1.0) are Galactic star forming regions in the southern 
sky, which are generally less-studied. Radio measurements indicate that IRAS 10049-5657 
is an extended HII region (Hill, 1968; Manchester, 1969). 

A number of distance estimates to IRAS 10049-5657 can be found in the literature, rang- 
ing from 5.1 to 7.1 kpc. Here we use the distance of 6.3 kpc estimated by Caswell & Haynes 
(1987) based on radio recombination line measurements (for R© = 8.5 kpc and Galactic 
rotation velocity = 220 km s~^ at R©). The distance to IRAS 10031-5632 is estimated to be 
3.7 kpc (Caswell & Haynes, 1987). 

Recently, it has been concluded that IRAS 10049-5657 complex harbours a very massive 
OB star cluster (Bik et al. 2005; Hanson et al., 2003). Bik et al. (2005) have carried out 
near-infrared K-band spectroscopy of few members of this cluster and find two very massive 
(03-04) stars here. IRAS 10049-5657 has been studied as a part of surveys for the search 
of emission lines including masers. Formaldehyde absorption has been detected towards this 
source at 4.8 GHz (Whiteoak and Gardner, 1974) and 14.5 GHz (Gardner and Whiteoak, 
1984). Whiteoak et al (1982) detected CO (1-0) line emission from this source using the 
4-m radio telescope of CSIRO. Search for methanol transition (Peng and Whiteoak, 1992), 
methanol maser (Schutte et al, 1993) and OH maser (Cohen et al, 1995) close to this source 
have led to negative results. IRAS 10049-5657 has been imaged by Puchalla et al. (2002) 
at 42 GHz using the Mobile Anisotropy Telescope on Cerro Toco (MAT/TOCO) and an 
integrated fiux of 22.2 ±2.1 Jy within a 0.3° circular beam has been obtained by them. Both 
these regions, IRAS 10049-5657 and IRAS 10031-5632 have been studied as a part of the 
Parkes-MIT-NRAO (PMN) survey (Kuchar and Clark, 1997). They find the peak radio fiux 
density of IRAS 10049-5657 and IRAS 10031-5632 at 5 GHz to be 25.7 Jy/beam and 1.0 



^Presently at INAF-Osscrvatorio Astrofisico di Arcetri, Largo E. Fermi 5, 50125 Firenze, Italy 
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Jy/beam, respectively, where the beam is ~ 4'. 9. CS (2-1) hne emission has been observed 
close to both these regions (Bronfman et al., 1996) using the SEST telescope. The IRAS-LRS 
(Low Resolution Spectrometer) spectrum of IRAS 10031-5632 shows strong [Ne II] emission 
line at 12.8 /xm (de Muizon et al., 1990; Simpson & Rubin, 1990), a relatively weak [S III] 
emission line at 18.7 /im as well as emission in the Unidentified Infrared Bands (UIBs) at 
7.7, 8.6 and 11.3 /im (Zavagno et al., 1992; de Muizon et al., 1990). Water vapor maser 
(Braz et al. 1989) as well as methanol maser (Schutte et al. 1993) have not been found close 
to IRAS 10031-5632. 

The less known southern Galactic massive star forming regions are being studied under 
a long-term program which involves observing these sources in far infrared (Vig et al. 2007, 
Ojha et al. 2002, Karnik et al. 2001, Ghosh et al. 2000, Verma et al. 1994). In this paper, 
we present a systematic study of the star forming regions associated with IRAS 10049-5657 
and IRAS 10031-5632. The star forming region associated with IRAS 10049-5657 is believed 
to harbour a cluster of very massive stars. IRAS 10031-5632 is, on the other hand, a young 
star forming region with few members belonging to the cluster. We have carried out an 
infrared study of these southern Galactic HII regions in detail with the aim of understanding 
the energetics, physical sizes, the spatial distribution of interstellar dust and its temperature 
as well as the associated young clusters. Section 2 describes observations and other data sets 
used. Section 3 describes the results and in Section 4, radiative transfer modelling of these 
sources is presented. A comprehensive discussion of these sources is carried out in Section 5 
and a brief summary is presented in Section 6. 



2. Observations and data reduction 

2.1. Far infrared observations 

The Galactic star forming regions associated with IRAS 10049-5657 and IRAS 10031-5632 
have been observed using the two-band far infrared (FIR) photometer system at the Cassegrain 
focus of the TIFR 100 cm (f/8) balloon borne telescope. The observations were carried out 
during the balloon fiight from the TIFR Balloon Facility, Hyderabad in India (latitude 17°. 47 
north, longitude 78°. 57 east) on Feb 20, 1994. Details of the telescope and the observational 
procedure are given by Ghosh et al (1988). The two FIR bands use a pair of 2x3 composite 
silicon bolometer arrays, cooled to 0.3K by liquid ^He which view identical parts of the sky 
simultaneously. The field of view of each detector is 1'.6. The absolute positions were es- 
tablished from the detections of catalogued stars with an optical photometer located at the 
focal plane of the telescope, simultaneously during the FIR observations. The planet Jupiter 
was observed at the beginning as well as at the end of the fiight. The observations of Jupiter 
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were used for the absolute flux calibration of the two band FIR photometer (12 channels) 
as well as for the determination of the Point Spread Function (PSF). The spectral response 
of each band of the FIR photometer was determined in the laboratory using a Michelson 
interferometer and a Golay cell as a comparison detector. The two FIR wavebands will be 
referred to as 150 and 210 /im bands corresponding to the Ag// for 30 K source with A~^ 
emissivity law. 

The region around IRAS 10049-5657 (~ 32' x 20') and IRAS 10031-5632 (~ 25' x 18') 
was mapped by raster scanning the region of the sky in cross-elevation with steps in elevation 
at the end of each scan. The FIR signals were gridded into a matrix with a pixel size of 
0'.3 X 0'.3. The deconvolution of the observed chopped signal matrix was carried out using 
the Maximum Entropy Method similar to that of Gull & Daniell (1978) (for details see Ghosh 
et al, 1988). An angular resolution of ~ 1' has been achieved in the FIR maps using this 
method. 



2.2. Other datasets used 

2.2.1. IRAS 

The data from the Infrared Astronomical Satellite (IRAS) survey in the four bands 
(12, 25, 60 and 100 /im) for IRAS 10049-5657 and IRAS 10031-5632 were HIRES (High 
Resolution processing using Maximum Correlation Method; Aumann et al., 1990) processed 
at the Infrared Processing and Analysis Center (IPAC, Caltech) to obtain high angular 
resolution maps. The flux densities of the sources within a circular region of diameter 3' 
(centered on the peak) have been extracted from these images. IRAS 10031-5632 appears in 
IRAS-LRS Catalog (IRAS Science Team, 1986) while the LRS spectrum of IRAS 10049-5657 
is presented by Volk & Cohen (1989). These spectra, in the wavelength range 8-22 /im, and 
the flux densities have been used for constructing the spectral energy distributions (SEDs). 
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2.2.2. MSX 



The Midcourse Space Experimenlo (MSX) surveyed the entire Galactic plane within 
|b| < 5° in four mid infrared wavebands: 8.3, 12.1, 14.7 and 21.3 /zm at a spatial resolution 
of ~ 18". 3 (Price et al., 2001). The panoramic images of the Galactic plane survey of MSX 
were taken from IPAC. The images of IRAS 10049-5657 and IRAS 10031-5632 were used to 
extract the sources and obtain the flux densities within a circular region of diameter 3' in 
order to construct the SEDs. Point sources close to these star forming regions have been 
selected from MSX Point Source Catalog Version 2.3 (Egan et al., 2003) and cross-correlated 
with 2MASS sources (see Sect. 2.2.3). 



2.2.3. 2MASS 

The point sources around the regions IRAS 10049-5657 and IRAS 10031-5632 were se- 
lected from the Two Micron All Sky Surve}a (2MASS) Point Source Catalog (PSC). The 
2MASS PSC is complete down to J < 15.8, H < 15.1 and K, < 14.3 mag for S/N> 10, 
in the absence of confusion. The 2MASS sources used in this study are those with good 
photometric quality (rdfig = 1 — 3). The J, H and K^ magnitudes of the selected sources 
have been used to construct colour-magnitude (CM) and colour-colour (CC) diagrams which 
have been used to study the embedded clusters in these regions. The JHK^ magnitudes and 
images were taken from IPAC. 



2.2.4. SUMSS 

The Sydney University Molonglo Sky Survey (SUMSS) is a radio imaging survey of the 
sky south of declination —30° (Bock et al., 1999). This survey uses the Molonglo Observatory 



^ This research made use of data products from the Midcourse Space Experiment. Processing of the data 
was funded by the Balhstic Missile Defense Organization with additional support from NASA Office of Space 
Science. This research has also made use of the NASA/ IPAC Infrared Science Archive, which is operated 
by the Jet Propulsion Laboratory, Caltech, under contract with the NASA. 

■^This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project 
of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of 
Technology, funded by the NASA and the NSF. 
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Synthesis Telescopcl (MOST), operating at 843 MHz with a bandwidth of 3 MHz. The 
synthesised beam size is ~ 43" x 43"cosec|(5|. The SUMSS radio images of IRAS 10049-5657 
and IRAS 10031-5632 have been used to study the distribution of ionised gas around this 
region. 



3. Results 

In this section, we present the results obtained from the observations as well as from the 
available data. As our study includes a wide wavelength range (near infrared to radio), we 
classify our results according to those for the dust and gas components (interstellar medium) 
versus the stellar component for both the regions, IRAS 10049-5657 and IRAS 10031-5632. 



3.1. IRAS 10049-5657 

3.1.1. Interstellar medium 

The deconvolved maps of IRAS 10049-5657 at 150 and 210 /^m are presented in Fig. [H 
Figured] (left) shows emission from the IRAS 10049-5657 region at 150 /im while Fig. [1] (right) 
shows the emission at 210 /im, from the complete region (32' x 20') mapped by the telescope. 
The angular size (50% contour level) of IRAS 10049-5657 at 150 and 210 /xm is ~ 2'. The 
FIR emission in these maps samples the cold dust around these regions. IRAS 10049-5657 
is well resolved in both the bands. The FIR emission in this region shows extension towards 
north-east in both the bands. This emission towards the north-east extends to a larger scale 
in the 210 /xm map in Fig. [T] (right). Taking advantage of the simultaneous observations 
in the two bands, with almost identical field of view, we have generated maps of the dust 
temperature (Td) and optical depth at 200 /im (r2oo)- The dust temperature T(150/210) 
and optical depth (r2oo) maps are shown in Fig. [2l For these maps, we have assumed dust 
emissivity of the form ex oc A~^. The temperature distribution shows plateaus of maximum 
(60 K) towards north-west and south-east of the FIR peak. We have detected dust as cold 
as 20 K. The peak optical depth at 200 /im is determined to be 0.007 and is located at the 
position of the peak of FIR emission. Another peak is seen to the west of the maximum 
optical depth with an extension corresponding to north-east. The IRAS-HIRES maps of the 
region around IRAS 10049-5657 at 12, 25, 60 and 100 /xm are shown in Fig. [31 The achieved 



''The MOST is operated by the University of Sydney and supported in parts by grants from the Austrahan 
Research CounciL 
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angular resolutions are 0'.7 x 0'.5 at 12 and 25 yum, 1'.4 x I'.l at 60 fim, 2M x 1'.9 at 100 /xm. 
The flux densities of IRAS 10049-5657 integrated over 3' diameter, centered on the peak, 
from the TIFR and IRAS-HIRES maps are listed in Table [TJ 

We have modelled thermal continuum emission from interstellar dust alongwith the 
emission in Unidentified Infrared Bands (UIBs), using mid-IR data from the MSX Galactic 
Plane Survey in the 8.3, 12.1, 14.7 and 21.3 /im bands. This has been carried out using 
the scheme developed by Ghosh & Ojha (2002). In this scheme, the emission from each 
(6" X 6") pixel in the MSX images is modelled to be a combination of two components: (i) 
thermal continuum from the warm dust grains (gray body) and (ii) the emission from the 
UIB features falling within the MSX band. The scheme assumes that dust emissivity follows 
the power law of the form ex oc A^^ and the total radiance due to UIBs in the 12 /xm band is 
proportional to that in the 8 /xm band. The spatial distribution of UIB emission predicted 
by this scheme is presented in Fig. HI The morphology of UIB emission shows extensions 
towards the south-west near the peak as well as due north-west at the fainter levels. The 
peak strength of the modelled UIB emission is 8.9 x 10^^ W m~^ Sr~^ and is close (~ 16") 
to the IRAS position. Table [D also lists the flux densities of IRAS 10049-5657 from MSX 
maps integrated within a circular region of diameter 3'. 

The SUMSS radio continuum emission from the region around IRAS 10049-5657 at 843 
MHz is shown in Fig. [51 A dynamic range of ~ 700 is achieved in this region (peak flux 
is 4.3 Jy/beam; rms noise is ~ 6 mJy/beam). The radio emission peaks at (0:2000 = iC^ 
06™ 41.15^ 62000 = -57° 12' 40.0"). The integrated radio flux density is ~ 42.4 Jy over 43.6 
arcmin^ . 



3.1.2. Stellar component 

The distribution of near-infrared sources (selected from 2MASS PSC) in the region 
around IRAS 10049-5657 was investigated. A higher density of sources close to the IRAS 
10049-5657 region compared to neighbouring regions implied an embedded cluster. We have 
used 2MASS sources to study the nature of embedded cluster. We flrst estimate the cluster 
radius. For this, we select a large region of radius ~ 300" around the IRAS peak. To account 
for the contribution from the fleld stars (for background determination), we select a control 
fleld which is 20' to the east of IRAS 10049-5657. The centre of the cluster was estimated by 
convolving a gaussian with the stellar distribution and taking the point of maximum density 
as the centre. To determine the radial profile, the cluster region was divided into a number 
of concentric annuli with respect to the cluster centre. The surface number density of stars 
were obtained by counting them in each 7" annulus and dividing by the annulus area. The 



King's model, f{r), and the inverse radius model, g{r), of the following functional forms are 
fitted to the surface density radial profile. 

/(r) = Ofc + 



1 + (r/rcky 



g[r) = ai + 

r 

Here, a^ and a^ are the fitted background constants for the King's and inverse radius profile, 

respectively. For the King's profile, fo represents the core concentration at radius zero and r^k 

is the core radius, oq represents the stellar density of the inverse profile at the core. The radial 

profile of the observed star density as well as the fits are shown in Fig. El A clear gradient 

in surface number density distribution confirms the existence of clustering (crowdiness) at 

centre of the region. The radial profile of the cluster merges with the background field at 

~ 65", yielding the extent of the cluster to be ~ 2 pc. The background level as estimated 

from the control field is ~3 stars pc^^, which agrees well with the values of 3.3 ± 0.6 and 

2.6 ±0.2 stars pc~^ yielded by the inverse radius and King's profile fitting, respectively. Note 

that these estimates of the stellar density towards the HII region when compared with the 

adjacent fields are unlikely to be affected by extinction (i.e. background stars) since the 

cluster is fairly distant. 

Within a cluster radius of 65" around the cluster centre (^2000 = 10^ 06'" 38.19"*, (^2000 
= -57° 12' 25.2"), 89 sources were found. Of these, 45 were detected in all the three JHKs 
bands of 2MASS. Forty- four sources have been detected in either H and K^ bands or only 
Kg band. We have discussed the nature of stellar populations in this region using the CM 
( J-H vs J) and the CC ( J-H vs H-K) diagrams which are shown in Fig. [71 We have assumed 
extinction values of Aj/Ay = 0.282, A^/Ay = 0.175 and AxJAy = 0.112 from Rieke & 
Lebofsky (1985). All the 2MASS magnitudes as well as the curves are in the Bessel & Brett 
(1988) system. In the CM diagram, the nearly vertical solid lines from left to right (with 
increasing J-H) represent the ZAMS curves (for a distance of 6.3 kpc) reddened by Ay = 0, 
15 and 30 mag, respectively. The slanting lines joining them trace the reddening vectors of 
these ZAMS stars. In the CC diagram, the locii of the main sequence and giant branches 
are shown by the solid and dotted lines, respectively. The short-dash line represents the 
locus of classical T-Tauri stars (Meyer et al., 1997). The parallel dot-dash straight lines 
follow their reddening vectors. The long-dash line represents the locus of Herbig Ae/Be 
stars (Lada and Adams, 1992). We classify the CC diagram into three regions, as shown 
in the figure (see Ojha et al., 2004). The 'F' region is considered to be the region where 
field stars (main sequence stars, giants). Class HI and Class II objects with small infrared 
excess are located. The 'T' sources are located redward of region 'F' but blueward of the 
reddening line projected from the red end of the T-Tauri locus of Meyer et al. (1997). These 
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sources may be considered to be mostly Class II with large NIR excess (Classical T Tauri 
like-stars) and/or extremely reddened early type ZAMS stars, having excess emission in K- 
band (e.g., Blum, Damineli & Conti, 2001 and references therein). The 'P' region is where 
the protostar-like Class I objects and Herbig Ae/Be stars are mostly located. In Fig. [7] (left) 
CM diagram, the 10 sources lying above the ZAMS curve of spectral type 09 are shown as 
asterisk symbols while 13 sources with infrared excess are shown as open circles. It is likely 
that few of these objects are foreground objects or bright background giants. The source 
lying above the ZAMS curve of 09 and lying in the 'T' region is shown as a solid triangle. 
The remaining sources are represented by plus-symbols. It is important to note that the 
CM and CC diagrams are useful tools for estimating the approximate nature of the stellar 
populations within the cluster in a statistical sense. 



3.2. IRAS 10031-5632 

3.2.1. Interstellar medium 

The intensity maps of the region around IRAS 10031-5632 in the 150 and 210 /xm FIR 
bands are shown in Fig. [HI The emission at 210 ^m. from the entire mapped region (25' x 18') 
is shown in Fig. [8] (right). The isophotes in both the maps have been displayed upto 1% level 
of the peak intensities. A detailed investigation of these maps reveals that IRAS 10031-5632 
is unresolved at 150 jum but is extended at 210 /xm map. At fainter levels, we observe 
large-scale extended emission, particularly towards the south-east. The HIRES-processed 
IRAS maps for the corresponding region around IRAS 10031-5632 are shown in Fig. O The 
angular resolution achieved in these maps are 0'.9 x 0'.5 at 12 and 25 /im, 1'.4 x I'.l at 60 
/xm, 2'.1 X 1'.9 at 100 ^m.. The peak position and flux density details of IRAS 10031-5632 
are given in Table [T] The emission in the mid infrared bands of MSX has been used to model 
the the peak UIB emission which is 6.5 x 10~^ W m~^ Sr~^. 

The SUMSS radio emission in region around IRAS 10031-5632 at 843 MHz is shown in 
Fig. [TOl The dynamic range of the radio map is ~ 300 (peak flux is 652 mJy/beam; rms 
noise is ~ 2 mJy/beam). The peak of radio emission is at 02000 = lO'^ 04™ 56.10"*, ^2000 = 
-56° 46' 38.0". The integrated flux density is 1.1 Jy over 4.9 arcmin^. 



3.2.2. Stellar component 

We have used the 2MASS sources in the vicinity of IRAS 10031-5632 to study the stellar 
populations here. In a circular region of radius 1' around the IRAS peak, 27 sources have 
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been detected in JHK^ bands while 14 sources are detected in either H and K^ bands or 
only K^ band. The CM (H-K vs K) and CC (J-H vs H-K) diagrams of the sources detected 
in all the three bands are shown in Fig. [TTl The ZAMS curves in the CM diagram are 
for a distance of 3.7 kpc reddened by Ay = 0, 15 and 30 mag, respectively. As in the 
case of IRAS 10049-5657, all the 2MASS magnitudes as well as the curves are in the Bessel 
& Brett (1988) system. The source lying above the ZAMS curve of spectral type 09 is 
represented by an asterisk while the filled triangle represents a source lying above the ZAMS 
curve of spectral type 09 and having an infrared excess. The 7 sources detected only in 
HKs bands are shown as open squares in the figure. The remaining sources are denoted by 
crosses. An investigation into the distribution of 2MASS PSC sources in the region around 
IRAS 10031-5632 shows a few sources grouped together near the IRAS peak. For this cluster, 
we have been unable to determine the cluster radius owing to the low surface density of stars 
(very few excess stars close to the IRAS peak). 



4. Radiative Transfer Modelling 

In an attempt to obtain a self-consistent picture of these star forming regions using 
all the available data as well as to extract important physical parameters, we have carried 
out radiative transfer modelling of IRAS 10049-5657 and IRAS 10031-5632. The symmetric 
morphology of contours near the centrally located peak in the FIR maps supports the one- 
dimensional treatment of the radiative transfer modelling. 



4.1. Continuum emission 

Each star forming region is modelled as a spherically symmetric cloud of gas (hydrogen) 
and dust, powered by a centrally embedded source (single or cluster of ZAMS stars). The 
cloud is immersed in an average interstellar radiation field. The radiative transfer equations 
have been solved assuming a two-point boundary condition for the spherical cloud. The gas 
exists throughout the cloud, i.e. from the stellar surface to the edge of the cloud. The dust 
exists in a spherical shell with a cavity at the centre. This is because close to the exciting 
source(s), the dust grains are destroyed when exposed to excessive radiative heating. Two 
commonly used types of interstellar dust are explored; the DL type (Draine & Lee, 1984) 
and the MMP type (Mezger, Mathis & Panagia, 1982). In the spherical shell where gas 
and dust co-exist, the gas-to-dust ratio is held constant. Other details of this self-consistent 
scheme are given in Mookerjea & Ghosh (1999). The parameters explored are as follows : 
the geometrical dimensions of the cloud (primarily the radius of the dust-free cavity; the 
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outer diameter is guided by the observed angular extent and the distance), radial density 
distribution law (three power law exponents : n(r) ~ r*^, r~^, or r~^), radial optical depth 
due to dust, the nature of embedded source/ (s) (single ZAMS star or a cluster, consistent 
with the total observed luminosity), relative abundances of different grain types (silicate, 
graphite), and the gas-to-dust ratio. The observational constraints include the SED due 
to thermal emission from the dust component, angular sizes at different wavelengths and 
radio continuum emission from the H II region. The parameters corresponding to the best 
fit models for the two sources are presented in Table 2. While no claim is made about the 
uniqueness of these parameter-sets, the following comments in support of the robustness of 
the results are in order : - a constant density radial profile provides the best fit for both 
regions and the r~^ or r~^ density profiles are conclusively ruled out by the observed SED; 
- inner (dust-free) cavity radii lower than those derived result in short wavelength fluxes 
significantly higher than those observed and are thus ruled out; - radio continuum emission 
is quite sensitive to the nature of the central energy source, viz., single star versus a cluster. 



4.2. Line emission from gas 

We have also modelled the infrared nebular /ionic fine structure line emission from gas 
around IRAS 10049-5657 and IRAS 10031-5632, using a sophisticated scheme which uses 
more details of the interstellar gas component. This scheme developed by Mookerjea & 
Ghosh (1999), uses the photoionisation code CLOUDY (Ferland, 1996), which solves for sta- 
tistical and thermal equilibrium by balancing various ionization-neutralization processes as 
well as heating-cooling processes. The scheme additionally includes : (a) the exact structure 
of the cloud / H II region (viz., central dust-free cavity) and (b) absorption effects of the 
dust component on the emergent line intensities. Typical H II region abundance of the gas 
component has been taken into consideration. The overall structure of the cloud is defined 
by the parameters that provide the best fit to the continuum spectral energy distribution 
(see subsection 4.1). The first part of the calculations involve the pure gas inner shell and 
the emerging spectrum comprises of continuum as well as line emission. This emergent con- 
tinuum from the inner shell provides the inner surface boundary condition for the second 
shell comprising of gas and dust. The line emission from the inner shell is transported out- 
wards through extinction by the dust column in the second shell. Similarly, line emission 
originating from the second shell is transported considering the absorption effects of the dust 
grains lying between the emission zone and the outer surface of the second shell. The finally 
predicted emergent line luminosity includes both these components. A total of 27 spectral 
lines in the wavelength range 2.5 — 200 /xm have been considered. To compare the predictions 
of the model with the spectral lines detected in the IRAS-LRS spectra of IRAS 10049-5657 
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and IRAS 10031-5632, we convolve the model predicted spectrum with IRAS-LRS instru- 
ment profile in the wavelength range 8-22 yum. However, to show the complete spectrum 
for the wavelength range 2.5-200 /im, the spectral lines have been convolved with ISO-SWS 
and LWS typical spectral resolutions of 1000 for 2.5 < A < 12 /xm, 20000 for 12 < A < 45 
/im, 8100 for 45 < A < 75 /im and 6800 for 75 < A < 200 /xm while predicting the expected 
emergent spectrum. 



4.3. IRAS 10049-5657 

The SED for IRAS 10049-5657 is constructed using the flux densities in the two TIFR 
bands, the four IRAS bands (from HIRES maps), the IRAS-LRS data and the four MSX 
bands. The flux densities used in the SED are the fluxes integrated over a circular region of 
diameter 3' around the peak. The total luminosity is 1 x 10^ L© for a distance of 6.3 kpc. 
Fig. [12] (left) shows the observed SED and the predicted spectrum from the best fit model. 
This best fit model implies a uniform dust and gas density distribution with the embedded 
energy source as a single ZAMS star of type 05-04 (see Table 2). The inner cloud dust 
radius is 0.008 pc while the outer cloud radius is 4.2 pc. The optical depth at 100 fim. is 
0.002. The radius of the ionised gas from the model is 2.8 pc. The observed angular sizes are 
explained by this model. The measured radio flux density, ~ 25 Jy at 843 MHz is obtained 
by integrating within a circular region of 2.8 pc around the radio peak. This flux is quite 
high and cannot be explained by the model as the predicted radio flux from the model is 4 
Jy for gas-to-dust ratio of 100:1 by mass. The dust composition for Si:Gr is 11:89 for the 
DL type of dust. The total dust mass is determined to be 12 Mq. 

The modelling of line emission for IRAS 10049-5657 has been carried out using the 
physical sizes obtained from the modelling of continuum emission. In all, 21 nebular/ionic 
lines satisfy our detectability criterion (power in the line be at least 1% of the power in the 
neighbouring continuum). The wavelengths and luminosities of these lines are presented in 
Table [31 The modelled spectrum, including the lines from the 10 elements considered as well 
as the continuum predicted by this model is shown in Fig. [12] (right). 

We have looked for signatures of the ionic lines generally found in HII regions in the 
IRAS-LRS spectrum using a simple method. This method assumes the line signal to be 
present within data points, I^nne, corresponding to one resolution element of LRS, centered 
at the expected line position. The local spectral baseline is estimated by interpolating 
(power law) data points in one resolution element each on both the red & blue sides of 
Diine- The estimated baseline is then subtracted from Dune to obtain the line emission. For 
IRAS 10049-5657, the 'brighter' lines: [Ne II] at 12.8 /xm, [Ne HI] at 15.5 /im, [S HI] at 18.7 
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/iin and [Ar III] at 21.8 fun are detected. The luminosities of these hnes as well as the ratios 
of line-to-continuum are listed in Table [31 

We compare the luminosities as well as the ratios of line-to-continuum of the predicted 
lines from the model (convolved with IRAS-LRS instrument profile) with respect to those 
observed. The line-to-continuum ratios for the 'LRS-convolved' model are listed in Table [31 
We find that the ratios of luminosities (model/observations) for the lines [Ne II], [Ne III] 
and [S III] agree within a factor of 4. For the [Ar III] line, we find that the observed value 
is ~ 13 times larger than the modelled value. However, it is to be noted that [Ar III] line is 
close to the edge of the wavelength range of the spectrometer and could have instrumental 
uncertainities. The hnes [Ar III] at 9.0 /xm and [S IV] at 10.5 /xm are barely detected although 
the model predicts them to be bright. 



4.4. IRAS 10031-5632 

The SED for IRAS 10031-5632 has been constructed using flux densities (integrated 
over a circular region of 3' diameter centred on the peak) from the MSX, IRAS-HIRES as 
well as the TIER maps. The IRAS-LRS spectrum has also been used in the construction 
of the SED. The spectrum shows a strong silicate feature. Eor IRAS 10031-5632, the total 
luminosity is 4.2 x 10^ Lq (d~ 3.7 kpc). We have used a single ZAMS star of spectral type 
09 as the centrally exciting source. The best fit radiative transfer model along with the 
observed SED is shown in Fig. [131 (left), and the parameters of the best fit model are given 
in Table [21 This model implies a uniform density distribution of gas and dust. The optical 
depth at 100 /xm (from the model) is 0.07. The ratio of silicates and graphite dust grains is 
estimated to be 62:38. The radio flux predicted by the model at 843 MHz is 0.07 Jy for a 
gas-to-dust ratio of 100 by mass and radius of ionised region is 0.03 pc. This is comparable 
to the measured radio flux density of 0.04 Jy at 843 MHz (obtained by integrating the flux 
within a circular region of radius 0.03 pc centred on the radio peak). A dust mass of 24 Mq 
is obtained from the model. 

For IRAS 10031-5632, 17 nebular /ionic lines satisfy the detectability criterion. The 
luminosities of the lines as well as the ratio of luminosities of each line with respect to 
the continuum (for the model convolved with the LRS spectral resolution) are also listed 
in Table [H Figure [T31 (right) shows the emerging spectrum from the model. In addition, 
following the method in Sect. 4.3, we have searched for lines in the LRS spectrum of 
IRAS 10031-5632. The [Ne II] line at 12.8 yum is clearly detected with a luminosity of 
30 Lq. Among the other lines in the range of LRS, we find that [Ar HI] at 9.0 fim and 
[S III] at 18.7 fim also show detections. A comparison of the luminosites of the detected 
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lines with the predicted values shows that the observed luminosities are much higher than 
the modelled ones. Table H] also lists the ratio of line-to-continuum for these lines. A 
comparison of these ratios shows that the model ratios are a factor of ~ 100 lower than 
those observed. This can be attributed to the elevated level of the continuum from the 
model. A comparison of the model SED with the observed one in Fig. [13] (left) shows that the 
radiative transfer model in the mid infrared overestimates the observed LRS spectrum. Much 
larger discrepancies between the observed and predicted line emissions for IRAS 10031-5632, 
may be understood as follows. Our scheme for predicting the line emission utilizes the 
description of the cloud, which is obtained from modelling the continuum emission from 
the dust component distributed throughout the cloud. In contrast, emission of the specific 
ionic lines under discussion, viz., [S III], [Ne II] and [Ar III] are expected to originate from 
the innermost part of the H II region (due to high ionization potentials). In addition, the 
lower excitation of the central source (ZAMS 09) for this source, makes these line emissions 
more sensitive to the precise physical details (e.g. inhomogeneities like dumpiness) in the 
immediate vicinity of the star. Since the emission of forbidden lines depends on square of 
the local density, denser clumps in an inhomogeneous medium will show enhanced emission 
compared to an equivalent uniform medium. Hence, a clumpy medium around the exciting 
star in IRAS 10031-5632, could help explain the observed higher nebular line luminosities. 
This scenario is also consistent with the higher radio continuum emission observed. 



5. Discussion 

5.1. IRAS 10049-5657 

Using the FIR map of IRAS 10049-5657 (Fig.[T]), we see an extension of the dust emission 
towards north-east. This is also seen in the 60 and 100 /xm emission IRAS-HIRES maps of 
this region, shown in Fig. [31 The IRAS-HIRES mid infrared emission from warm dust, 
however, shows an extension towards the west. This western extension is also clearly seen 
in the MSX maps. The flux densities from the TIER maps at 150 and 210 /im have been 
used to compute the mass of dust using the formulation of Hildebrand (1983) and Sandell 
(2000). For a temperature of 30 K, the dust mass obtained is ~ 21 Mq. As can be seen from 
Fig. [S], a north-east extension seen also in the cold dust emission is observed for the ionised 
gas also. The total flux of 42.4 Jy at 843 MHz can be compared with the flux density of 40 
Jy (beam ~ 14') at 1410 MHz (Manchester, 1969). 

From the CM (J-H vs J) diagram of the 2MASS sources (Fig. [7] left) within the cluster 
radius (~ 65"), we see that there are 11 sources lying above the ZAMS curve of 09. These 11 
sources are designated as IRA-1, ... , IRA-11 and a list of their positions and magnitudes is 
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given in Table [51 However, it is important to note that the spectral types inferred from CM 
diagram are the earliest possible spectral type (upper limits) when the stars have infrared 
excess. Sources with infrared excess can be found from the CC (J-H vs H-K) diagram. In 
the CC diagram (see Fig. [7]), 8 'asterisks' lie in the band occupied by reddened ZAMS stars. 
One source, IRA-7 (shown as a triangle lying in the 'T' region), shows an infrared excess 
while IRA-10 lies to the left of the reddening band of the ZAMS objects (drawn from the top 
of the main sequence branch) in the CC diagram. The upper left part of the CC diagram is 
not an allowed region for young stellar objects (Lada & Adams, 1992). IRA-10 is faint in J 
band (16.9 ±0.2 mag) with larger errors though it is relatively brighter in the H (14.2 ±0.02 
mag) and Kj, (13.2 ± 0.04) bands as compared to the J band. One possibility to explain 
the near infrared colours of IRA-10 as well as its position in the CC diagram is that IRA- 
10 could comprise two or more unresolved sources. The sources lying near the unreddened 
main sequence in the CC diagram are possibly foreground sources not associated with the 
cluster. A small but significant fraction of cluster stars lie outside and to the right of the 
reddening band of the ZAMS stars. These are mostly young stellar objects (YSOs) with 
intrinsic colour excess. Fourteen objects lie in this infrared excess zone, i.e. in the 'T' and 
'P') regions. By dereddening the stars (on the CC diagram) that fall within the reddening 
vectors encompassing the main sequence stars, we find the visual extinction (Ay) towards 
each star. The individual extinction values range from Ay ~ — 17 magnitudes. From the 
CM diagram, we find that the extinction values of most 'asterisks' lie at about Ay ~ 15 — 18 
magnitudes. 

Figure [H] shows the 2MASS K^-band image of the region around IRAS 10049-5647 in 
grayscale. The 2MASS image of this region shows diffuse emission apart from the sources 
(earlier than spectral type 09) clustered together. It is interesting to note the clustering of 
these sources close to the IRAS peak and the distribution of the other such stars in the NE- 
SW direction which is the direction of extension of the ionised gas as well as the cold dust. 
The morphology of emission in the UIB as well as the emission from warm dust in MSX 
bands also shows extension along the south-west direction. The brightest infrared source 
among the selected 2MASS sources is IRA-11. IRA-7 is closest (~ 2.7") to the IRAS peak, 
is of spectral type earlier than 05 and has an infrared excess. This is consistent with the 
spectral type determined by Bik et al. (2005) for IRA-7 (referred to as 10049nr411 in their 
paper) using high resolution K-band spectra. They find that it is of spectral type 03-04. 
Recent multi-epoch radial velocity measurements by Apai et al. (2007) show large amplitude 
variations in radial velocities pointing towards IRA-7 being a massive binary system (50 and 
20 Mq). Bik et al. (2005) have also carried out spectroscopy of IRA-3 (10049nr324) and 
find it to be of spectral type 03-04/05-06. This is consistent with its ZAMS spectral type, 
05-06, obtained from the CM diagram. It is interesting to note that the 'asterisk' closest 
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to the radio peak is IRA-8. The CM diagram indicates IRA-8 to be of ZAMS spectral type 
09-06. 

We have also cross-correlated the MSX PSC sources with those from 2MASS PSC which 
lie within the cluster radius. There are 2 such MSX PSC sources namely G282. 0341-01. 1810 
and G282.0176-01.1793 (hereafter Ml and M2, respectively; listed in Table [6]). Ml coincides 
(1.3") with an infrared excess source, ^mass-j 1006411^ -5712377 (J~ 16.0 ±0.2, H~ 14.9 ±0.2, 



Kg ~ 13.5±0.3). M2 coincides (2") with a source detected only in K^ band, 2mass-j 1006359 2- 
5711563 (K, ~ 14.1 ± 0.2 mag). From the 2MASS-MSX colour-colour diagram (Fai/Fg vs 
Fs/Fk] Lumsden et al. 2002, their Fig. 9), we find that these sources lie in the region 
generally covered by compact H II regions and massive young stellar objects. We can, 
therefore, conclude that these are young stars associated with the cluster. 



5.2. IRAS 10031-5632 

The dust emission around IRAS 10031-5632 is compact (unresolved) at 150 /xm, as well 
as at the IRAS-HIRES wavebands. IRAS 10031-5632 is not resolved in the SUMSS radio 
map also. However, IRAS 10031-5632 is barely resolved in the 210 /im map. Deconvolving 
the beam from 210 /xm image gives us ~ 0.9' as an estimate of its size. The dust mass 
obtained using the flux density at 210 yum from the TIFR map is ~ 7 Mq for a temperature 
of 30 K. However, from the radiative transfer modelling, we obtain a larger dust mass of 24 
Mq. The total radio flux density is 1.1 Jy. Considering that the radio flux is from a single 
ionising ZAMS star, we use the formulation of Schraml & Mezger (1969) and Panagia (1973) 
to estimate its spectral type. We flnd it to be of ZAMS spectral type 09-08.5. This can be 
compared with the ZAMS spectral type of 09 obtained from the FIR luminosity. 

Unlike the case of IRAS 10049-5657, we do not see a rich cluster with a large number of 
stars around IRAS 10031-5632. Rather a small group comprising of a few stars is seen close 
to this region. There are seven sources (labelled as IRB-1 to IRB-7) clustered close to the 
FIR (IRAS) peak; of which two, IRB-1 and IRB-5, are detected in all the three (J, H and 
Ks) bands. The 2MASS designation and flux details of these sources are given TableH From 
the CM diagram (H-K vs K; Fig. [IDleft), it is observed that IRB-1 is of ZAMS spectral type 
earlier than 05. It is to be noted that this spectral type of IRB-1 obtained from NIR study 
is the earlier possible spectral type (upper limit) as it shows an infrared excess. IRB-5 is 
however of later spectral type, B3-B2. It lies among a well-deflned group of stellar sources 
which have lower extinction values (Ay ~ 5 mag). From the CC diagram (J-H vs H-K), 
it is found to lie near the main-sequence curve. It is therefore possible that IRB-5 is not 
associated with the cluster but is a foreground source. The other flve sources are detected 
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in H and K^ or only in K^ bands. The source closest (~ 1.5") to the IRAS peak, IRB-3, is 
detected in H (15.1 ± 0.2 mag) and K^ (12.8 ± 0.1 mag) bands but not detected in the J 
band. From the CM diagram, it is found that IRB-3 is a heavily extincted early spectral 
type (09-06) ZAMS star. It is located at the peak of 843 MHz radio emission. The other 
four sources IRB-2, IRB-4, IRB-6 and IRB-7 are detected only in K^ band. This implies 
that these are deeply embedded objects. Thus, among the small number of cluster members, 
the majority of them are deeply embedded indicating that the star formation here is in an 
early stage. The 2MASS image of the region around IRAS 10031-5632 is shown in Fig. [T51 
The sources lying above the ZAMS curve of spectral type 09 are shown as asterisk symbols 
while the 'asterisk' having an infrared excess is shown as a solid triangle. The open squares 
represent sources detected either in H and K^ bands or only in K^ band. The other sources 
are represented with plus-symbols. The seven sources discussed above are labelled as "IRB" 
in the figure. It is interesting to note the presence of young objects (detected in H and K^ 
or only K^ bands) clustered near the IRAS peak. 

In the IRAS 10031-5632 region, we find that there is 1 MSX PSC source within 1' radius 
centered on the IRAS peak, G281. 5857-00.9706 (designated as M3 and listed in Table E]). 
M3 could possibly be associated with either IRB-3 (~ 2.8") or IRB-2 (~ 2.8"). In both the 
cases, we obtain the 2MASS-MSX colours (F21/F8 and F^/Fk, Lumsden et al. 2002) and 
observe that M3 lies in the general region covered by compact H II regions. 



6. Summary 

The massive star forming regions associated with IRAS 10049-5657 and IRAS 10031-5632 
have been studied using the infrared (near, mid and far) wavebands. The dust and gas envi- 
ronments as well as the stellar sources have been probed using data from TIFR balloon-borne 
telescope, MSX, SUMSS and 2MASS. The spatial distributions of far infrared emission from 
cold dust at 150 and 210 /im have been obtained alongwith maps of optical depth (r2oo) 
and colour temperature, T(150/210), Using MSX data, the emission from warm dust and 
UIBs in the IRAS 10049-5657 region has been studied. IRAS 10049-5657 region shows the 
presence of a rich cluster of OB stars which gives rise to strong radio continuum (~ 42 Jy 
at 843 MHz). The cluster radius is estimated to be ~ 65" and within the cluster radius, 
there are 11 2MASS sources lying above the ZAMS curve of spectral type 09, designated 
IRA-1 to IRA-11. One source, IRA-7 also has an infrared excess and is located closest to the 
IRAS peak. It is a massive star of early spectral type, earlier than 05 consistent with earlier 
studies. Unlike IRAS 10049-5657, the star forming region associated with IRAS 10031-5632, 
comprises of a very compact group of stars. The SUMSS radio flux at 843 MHz predicts 
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a radio spectral type of 09-08.5. The sources clustered near the IRAS peak comprise of 
young sources detected in H and K^ or only K^ bands. The 2MASS source detected in 
all the three bands, designated IRB-1, has an infrared excess and is found to be of ZAMS 
spectral type earlier than 05 from the 2MASS colour- magnitude diagram. Self-consistent 
radiative transfer modelling constrained by observations have been carried out for both these 
sources. The geometric details of the clouds, the dust-composition and optical depths, etc 
have been obtained from the best fit models. We have also carried out modelling of line 
emission from IRAS 10049-5657 and IRAS 10031-5632 using a scheme based on CLOUDY. 
The predictions of the model for emission in the ionic lines are closer to LRS detections for 
IRAS 10049-5657 than for IRAS 10031-5632. We speculate that for IRAS 10031-5632, some 
of the basic assumptions of the modelling scheme are not valid. 
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Table 1: Flux density details of IRAS 10049-5657 and IRAS 10031-5632 
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Table 2: Best-fit parameters of the radiative transfer models for IRAS 10049-5657 and 
IRAS 10031-5632. 
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Composition 
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0.03 
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10 
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DL 
DL 


11:89 
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12 
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" Luminosity of single ZAMS star used in the model 
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Table 3: Emergent line luminosities predicted by the model for IRAS 10049-5657. The line- 
to-continuum ratios for the model (convolved with the IRAS-LRS instrument profile) are 
also listed. Also included are the line luminosities and the ratio of line-to-continuum for the 
hues detected in the IRAS-LRS spectrum of IRAS 10049-5657. 
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Element 
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Ar III 
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263 0.04 


S III] 


18.69 
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Ne III 


15.57 
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- 
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- 
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- 
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- 
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Table 4: Emergent line luminosities predicted by the model for IRAS 10031-5632. The line- 
to-continuum ratios for the model (convolved with the IRAS-LRS instrument profile) are 
also listed. Also included are the line luminosities and the ratio of line-to-continuum for the 
lines detected in the IRAS-LRS spectrum of IRAS 10031-5632. 
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- 


S III 


33.50 


0.53 




- 


[Ar III] 


21.84 


0.10 


0.0002 


- 


S III 


18.69 


1.76 


0.0025 


16.6 0.07 


[Ne III 


15.57 


3.13 


0.004 


- 


[Ne II] 


12.82 


1.61 


0.005 


30 0.5 


[S IV] 


10.52 


0.27 


0.0006 


- 


[Ar III] 


8.99 


0.40 


0.0009 


9.3 0.12 


[Ar II] 


7.00 


0.17 




- 
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Table 5: Details of 2MASS PSC sources earlier than 09 around IRAS 10049-5657. 



2MASS PSC 


Name"' 


^2000 


(^2000 


J 


H 


K, 


designation 




(deg) 


(deg) 


(mag) 


(mag) 


(mag) 


J10063343-5713064 


IRA-1 


151.639294 


-57.218468 


13.18 ±0.05 


11.41 ±0.05 


10.66 ±0.04 


J10063408-5713079 


IRA-2 


151.642008 


-57.218887 


13.57 ±0.02 


12.28 ±0.02 


11.79 ±0.02 


J10063696-5712372 


IRA-3 


151.654013 


-57.210342 


14.46 ±0.03 


12.72 ±0.03 


11.84 ±0.02 


J10063849-5712287 


IRA-4 


151.660405 


-57.207985 


14.57 ±0.08 


12.90 ±0.12 


11.92 ±0.06 


J10063875-5712243 


IRA-5 


151.661484 


-57.206776 


14.56 ±0.09 


12.75 ±0.12 


11.80 ±0.10 


J10063940-5712198 


IRA-6 


151.664190 


-57.205513 


13.79 ±0.03 


11.97 ±0.05 


10.95 ±0.04 


J10063949-5712299 


IRA-7 


151.664564 


-57.208324 


12.89 ±0.04 


11.44 ±0.06 


10.36 ±0.04 


J10064088-5712327 


IRA-8 


151.670346 


-57.209087 


14.97 ±0.09 


13.20 ±0.08 


12.13 ±0.10 


J10064153-5712003 


IRA-9 


151.673054 


-57.200096 


14.77 ±0.03 


12.90 ±0.03 


11.89 ±0.04 


J1006421 1-5712041 


IRA- 10 


151.675482 


-57.201160 


16.88 ±0.18 


14.21 ±0.02 


13.17 ±0.04 


J10064217-5712231 


IRA- 11 


151.675725 


-57.206444 


11.33 ±0.02 


9.96 ±0.03 


9.38 ±0.03 



"" Short name used in the present work 



Table 6: Details of the MSX PSC sources close to IRAS 10049-5657 and IRAS 10031-5632 (see text for details). 



IRAS 


MSX PSC 


Name" 


tt2000 


^2000 


Fs 


Fl2 


Fu 


F21 


Association 


designation 




(deg) 


(deg) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


IRAS 10049-5657 


G282.0341-01.1810 


Ml 


151.6721 


-57.2105 


16.5 ±0.7 


56.9 ±2.8 


97.0 ±5.9 


195.1 ±11.7 


IRAS 10049-5657 


G282.0176-01.1793 


M2 


151.6493 


-57.1995 


3.2 ±0.1 


13.1 ±0.7 


23.5 ±1.4 


39.2 ±2.3 


IRAS 10031-5632 


G281.5857-00.9706 


M3 


151.2342 


-56.7765 


8.0 ±0.3 


14.7 ±0.7 


17.1 ±1.0 


93.2 ±5.6 



to 



Short name used in the present work 
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Table 7: Details of 2MASS PSC sources close to IRAS 10031-5632. 



2MASS PSC 


Name" 


^2000 


^2000 


J 


H 


K, 


designation 




(deg) 


(deg) 


(mag) 


(mag) 


(mag) 


J10045536-5646422 


lRB-1 


151.230694 


-56.778397 


15.25 ±0.09 


12.56 ±0.07 


10.70 ±0.04 


J10045591-5646340 


lRB-2 


151.232992 


-56.776119 


— 


— 


13.24 ±0.20 


J10045627-5646381 


lRB-3 


151.234494 


-56.777267 


— 


15.14 ±0.24 


12.83 ±0.13 


J10045634-5646226 


lRB-4 


151.234762 


-56.772972 


— 


— 


13.80 ±0.18 


J10045653-5646275 


lRB-5 


151.235575 


-56.774307 


12.14 ±0.03 


11.85 ±0.04 


11.69 ±0.06 


J10045658-5646324 


lRB-6 


151.235780 


-56.775673 


— 


— 


13.11 ±0.13 


J10045666-5646424 


lRB-7 


151.236113 


-56.778446 


— 


— 


14.02 ±0.11 



" Short name used in the present work 
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Fig. 1. — The intensity map for the region around IRAS 10049-5657 at 150 /im (left) and 
210 /im (right). While the 150 /xm map shows emission from IRAS 10049-5657 region, the 
210 /im map shows the complete region scanned by the far-infrared telescope. Contour levels 
are at 1, 5, 10, 20, 40, 60, 80, 95% of peak intensity of 624 Jy/sq arc min (left) and 292 
Jy/sq arc min (right). 
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Fig. 2. — The distribution of dust temperature T(150/210) (left), and optical depth at 200 
/im, T200, (right) from the region around IRAS 10049-5657 assuming a dust emissivity law 
of eA oc A-2. The isotherms of T(150/210) correspond to 20, 27, 30, 33, 40 and 60 K. The 
T200 contours represent 5, 10, 20, 30, 40, 50, 60, 70, 80 and 90 % of the peak value of 0.007. 
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Fig. 3. — The IRAS-HIRES intensity maps for the region covering IRAS 10049-5657 at 12 
lira (top left), 25 /xm (top right), 60 fim (bottom left) and 100 /xm (bottom right). The 
contours are at 1, 5, 10, 20, 40, 60, 80, 90 and 95 % of the peak value of 272 Jy arcmin^^, 
1053 Jy arcmin^^, 3367 Jy arcmin^^ and 1423 Jy arcmin^^ at 12, 25, 60 and 100 /xm, 
respectively. 
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Fig. 4. — The spatial distribution of total radiance in UIBs for the region around 
IRAS 10049-5657 modelled using the MSX images. The contour levels are at 5, 10, 20, 
30, 40, 50, 60, 65, 70, 80, 90, 95% of peak value of 8.9 x 10^^ W m^^ Sr~\ 
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Fig. 5. — The SUMSS radio flux density map for tlie region around IRAS 10049-5657 at 843 
MHz. The contour levels are at 0.1, 1, 5, 10, 20, 30, 40, 50, 60, 70, 80, 90, 95% of peak flux 
of 4.3 Jy/beam. The synthesized beam is 43" x 43". 
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Fig. 6. — The radial profile of the surface number density for the cluster associated with 
IRAS 10049-5657 in log-log scale. Also plotted are the two fitted models - King's model 
(solid) and inverse radius model (dotted). The horizontal, dashed line corresponds to the 
background field star level, which is ~3 stars pc~^. Statistical errors (corresponding to the 
Poisson noise) are shown. 
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Fig. 7. — Colour-magnitude (left) and colour-colour diagram (right) for sources detected in 
all the three 2MASS bands for the region around IRAS 10049-5657. In the colour-magnitude 
diagram, the nearly vertical solid lines from left to right represent the zero age main sequence 
(ZAMS) curves reddened by Ay = 0, 15 and 30 magnitudes, respectively. The slanting lines 
trace the reddening vectors of these ZAMS stars. In the colour-colour diagram, the locii of 
the main sequence and giant branches are shown by the solid and dotted line, respectively. 
The short-dash line respresents the locus of classical T Tauri stars. The two parallel dot-dash 
straight lines follow the reddening vectors drawn from the base of curves of giants, main- 
sequence stars and T- Tauri stars. The long-dash line represents the locus of Herbig Ae/Be 
stars. The colour-colour diagram is classified into three different regions, namely 'F', 'P' and 
'T' (see text for details). The asterisk symbols represent sources lying above the ZAMS curve 
of spectral type 09 from the colour-magnitude diagram. The open circles represent sources 
depicting an infrared excess in the colour-colour diagram. The solid triangle represents a 
source lying above the ZAMS curve of spectral type 09 as well as having an infrared excess 
i.e. lying in the 'T' region. The plus symbols represent the other sources. 



-33- 



-55" 50' L 




10 05 30^ 10" 04'" 50" 

Right Ascension (J2000) 




Right Ascension (J2000) 



Fig. 8. — The intensity map for the region around IRAS 10031-5632 at 150 /xm (left) and 
210 fira (right). While the 150 /xm map shows emission from IRAS 10049-5657 region, the 
210 /im map shows the complete region scanned by the far-infrared telescope. Contour levels 
are at 1, 5, 10, 20, 40, 60, 80, 95% of peak intensity of 712 Jy/sq arc min (left) and 220 
Jy/sq arc min (right). 
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Fig. 9.— The IRAS-HIRES intensity maps for the region covering IRAS 10031-5632 at 12 
lira (top left), 25 /xm (top right), 60 /im (bottom left) and 100 /xm (bottom right). The 
contours in the emission maps are at 1, 5, 20, 40, 60, 80, 90 and 95 % of the peak value of 
41.8 Jy arcmin"^, 774 Jy arcmin^^, 1370 Jy arcmin^^ and 486 Jy arcmin^^ at 12 and 25, 60 
and 100 /im respectively. 
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Fig. 10. — The SUMSS radio flux density map for tlie region around IRAS 10031-5632 at 
843 MHz. Tlie contour levels are at 0.5, 1, 5, 10, 20, 30, 40, 50, 60, 70, 80, 90, 95% of peak 
flux of 0.65 Jy/beam. The beam is 43" x 43". 
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Fig. 11. — Colour-magnitude (left) and colour-colour diagram (right) for sources detected in 
all the three 2MASS bands for the region around IRAS 10031-5632. In the colour-magnitude 
diagram, the nearly vertical solid lines from left to right represent the zero age main sequence 
(ZAMS) curves reddened by Ay = 0, 15 and 30 magnitudes, respectively. The slanting lines 
trace the reddening vectors of these ZAMS stars. In the colour-colour diagram, the locii of 
the main sequence and giant branches are shown by the solid and dotted line, respectively. 
The short-dash line respresents the locus of classical T Tauri stars. The three parallel dot- 
dash straight lines follow the reddening vectors drawn from the base of curves of giants, 
main-sequence stars and T- Tauri stars. The long-dash line represents the locus of Herbig 
Ae/Be stars. The colour-colour diagram is classified into three different regions, namely 
'F', 'P' and 'T' (see text for details). The asterisk symbols represent sources lying above 
the ZAMS curve of spectral type 09 from the colour-magnitude diagram. The open circles 
represent sources depicting an infrared excess in the colour-colour diagram. The solid triangle 
represents a source lying above the ZAMS curve of spectral type 09 as well as having an 
infrared excess. The plus symbols represent the other sources. The open squares in the CM 
diagram represent the 2MASS sources detected in only H and K^ bands. 
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Fig. 12. — (Left) Comparison of the spectral energy distribution from observations and the 
best fit radiative transfer model of IRAS 10049-5657 for continuum emission. The filled 
circles represent the TIFR and IRAS-HIRES data. The dashed lines are the IRAS-LRS 
spectra and the open squares represent the MSX data. The solid line presents the predicted 
emergent spectrum from our radiative transfer model. See text and Table [2] for details of 
model parameters. (Right) Emergent spectrum predicted by the scheme for line emission 
modelling for IRAS 10049-5657. The solid line represents the model convolved with the 
ISO spectral resolutions while the dot-dashed hne (shifted along Y-axis for better viewing) 
represents convolution with the LRS-instrument profile. 
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Fig. 13. — (Left) Comparison of the spectral energy distribution from observations and the 
best fit radiative transfer model of IRAS 10031-5632 for continuum emission. The filled 
circles represent the TIFR and IRAS-HIRES data. The dashed lines are the IRAS-LRS 
spectra and the open squares represent the MSX data. The solid line presents the predicted 
emergent spectrum from our radiative transfer model. See text and Table [2] for details of 
model parameters. (Right) Emergent spectrum predicted by the scheme for line emission 
modelling for IRAS 10031-5632. The solid line represents the model convolved with the 
ISO spectral resolutions while the dot-dashed hne (shifted along Y-axis for better viewing) 
represents convolution with the LRS-instrument profile. 
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Fig. 14. — The grayscale 2MASS K^-band image for the region around IRAS 10049-5657. 
The labelled axes are in J2000 coordinates. Overlaid on top are the SUMSS 843 MHz radio 
contours. The location of the IRAS peak is represented by a cross. Various 2MASS PSC 
sources are marked. The asterisk symbols represent sources lying above the ZAMS curve 
of spectral type 09. The open circles represent sources depicting an infrared excess in the 
colour-colour diagram. The filled triangles represent sources lying above the ZAMS curve of 
spectral type 09 as well as having an infrared excess. The other sources are marked by plus 
symbols. 
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Fig. 15. — The grayscale 2MASS K^-band image for the region around IRAS 10031-5632. 
The labelled axes are in J2000 coordinates. Overlaid on top are the SUMSS 843 MHz radio 
contours. The location of the IRAS peak is represented by a cross. Various 2MASS PSC 
sources are marked. The asterisk symbol represents the sources lying above the ZAMS curve 
of spectral type 09. The solid triangle represents the source lying above the ZAMS curve of 
spectral type 09 and having an infrared excess. The open squares represent sources detected 
in either H and Kg bands or only K^ band. The plus symbols represent the other sources. 



